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Abstract 



O" 1 In the framework of SUSYQM extended to deal with non-Hermitian Hamil- 

tonians, we analyze three sets of complex potentials with real spectra, recently 
derived by a potential algebraic approach based upon the complex Lie algebra 
j_j ■ sl(2,C). This extends to the complex domain the well-known relationship be- 

tween SUSYQM and potential algebras for Hermitian Hamiltonians, resulting 
from their common link with the factorization method and Darboux transfor- 
mations. In the same framework, we also generate for the first time a pair of 
elliptic partner potentials of Weierstrass p type, one of them being real and the 
other imaginary and PT symmetric. The latter turns out to be quasiexactly 
solvable with one known eigenvalue corresponding to a bound state. When 
the Weierstrass function degenerates to a hyperbolic one, the imaginary po- 
tential becomes PT non-symmetric and its known eigenvalue corresponds to 
an unbound state. 
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1 Introduction 



Recently there has been some interest 0, @, [5[ f§, [§, [7J] in studying PT-symmetric 
quantum mechanical systems. In quantum mechanics, the Hamiltonian of the under- 
lying system is usually assumed Hermitian ensuring a real energy spectrum. However 



ing the Schrodinger Hamiltonian to be invariant under the joint action of parity (P) 
and time reversal (T) transformations, one can still have a real spectrum of energy 
eigenvalues. Moreover, the overall normalizability of wave functions in many cases 
is not affected. In the literature, PT-symmetric schemes have been explored with 
respect to the complexification of several well-known potentials. Further, new ones 
have been searched for using a variety of techniques || f|, [5|, || [7]]. 

In this paper, we shall present results from a supersymmetric point of view. We 
shall show that the constraints furnished by the commutation relations of the sl(2,C) 
algebra admit of solutions that are consistent with supersymmetric intertwining 
relations. We shall also exploit the latter to obtain elliptic solutions of Weierstrass 
p type. Indeed we shall report here for the first time a PT-symmetric potential 
defined in terms of Weierstrass p function. 



2 N = 2 SUSYQM and Its Complexification Pro- 
cedure 



In order to set the notations, it would be useful to briefly recall the key features 
of the one-dimensional N = 2 SUSY quantum mechanics (SUSYQM). As is well 
known ||, the latter involves a pair of supercharges Q and Q\ related by Hermitian 
conjugation, and in terms of which the governing Hamiltonian H s is expressed as 



The supercharges Q and are fermionic in character and commute with H s : 



it has been conjectured JTJ that under less restrictive situations, namely by requir- 




(2.1) 



q 2 = (g f ) 2 = o, 



[Q,H S ] 



= Q\H a =0. 



(2.2) 
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A convenient way to deal with Q and is to adopt the representations 

Q = A®ct_, Q ] = A ] ®a +1 (2.3) 

where A is some linear differential operator and a± are combinations a± — <J\ ± icr 2 
of the Pauli matrices. A first-derivative realization of A and A\ namely 

A = ^- + W(x), ^ = -JL + w(x), (2.4) 
ctx ax 



yields the forms 



where is the so-called superpotential of the system. Note that fl2.5|) renders 

if, diagonal, 

ft = ( 1 ) • (2.6) 
As such we can factorize H± in the manner 

H + = A jf A = —^- + V M {x)-E, 
ax z 

H = AA^ = -JL + V ( --\x)-E, (2.7) 

at some arbitrary factorization energy E , with the partner potentials V^' related 
to W(x) through 

= W 2 tW + E. (2.8) 

It is easy to be convinced that the spectra of H + and H_ are alike except possibly 
for the ground state. In the exact SUSY case to which we shall restrict ourselves 
here, the ground state at vanishing energy is nondegenerate and is associated with 
H + or H , 

(x) = 0, 4 +) (x) = K exp (- J X W(t)dt\ , (2.9) 



or 



H_^ { ~\x) = 0, ipt } \x) = ^ ex P (f W(t)dtj , (2.10) 



according to whether J x W(t)dt — > +00 or —00 as x — > ±00. In ( |2.9| ) and ( |2.10| ), K 
is some normalization constant. 

The double degeneracy of the spectrum of H s can also be summarized by inter- 
twining H + and H_ according to 

AH + = H_A, H + A ] = A ] H_. (2.11) 

These relations follow from ( |2.7| ). 

To generate non-Hermitian potentials within SUSYQM [0, it is instructive to 
decompose the underlying superpotential W(x), the partner potentials V^(x), and 
the factorization energy E into a real and an imaginary part, namely 

W(x) = f(x)+ig(x), (2.12) 

V i+ \x) = vi t + \x) + iV I {+ \x), (2.13) 

V { -\x) = Vjr\x) +il / / ( " ) (x), (2.14) 

E = E R + iEj, (2.15) 

where /, g, V^, V^, Er, and Ej G R. All this leads to SUSY without Hermiticity: 
in particular, the supercharges are no longer related by Hermitian conjugation. As 
will be evident below, this presents no difficulty in so far as developing a theoretical 
framework is concerned. 

From (U) and (^T2|) -(^T5|), it follows that 



v£ +) = f 2 -g 2 -f' + E R , (2.16) 

V} +) = 2fg-g' + E I , (2.17) 

V { R - ] = f 2 -g 2 + f' + E R , (2.18) 

Vt ] = 2fg + g' + Ej. (2.19) 

These expressions are consistent with intertwining relationships. 

Since we will be interested only in a real energy spectrum, we can set Ej = 0. 
As such our basic relations correspond to 

V& +) = f 2 -g 2 -f' + E R , (2.20) 



V, (+) = 2fg-g', (2.21) 
V { r ] = f 2 -g 2 + f' + E R , (2.22) 
Vt } = 2fg + gf. (2.23) 

Observe that and Vr ^ are related by / — > —f, while Vj + ^ and Vj ^ are linked 
through g — > — y. In the following our task will be to analyze plausible solutions 
for the functions / and g pertaining to the set fl2.20|) - (|2.23|) , including the PT- 



symmetric ones. The main point to be noted is that by imposing an additional 
PT-symmetric restriction the Hermitian character of the intertwined Hamiltonians 
is lost. 

We shall first of all seek connections with the sl(2,C) potential algebraic approach 
to the construction of non-Hermitian Hamiltonians with real spectra. To this end 
we will show that our potentials ( |2.2(J )-(2.23) fit into such a scheme of complex po- 



tentials. In the next section, we therefore make a few remarks about the realization 
of the potential algebra sl(2,C). 

3 sl(2,C) Potential Algebra 

In Ref. ||, we made a detailed study of the sl(2,C) algebra. Its underlying commu- 
tation relations are 

[Jo, J±] = ±J±, [J+, J_] = -2J Q . (3.1) 
The generators Jq and J± can be given a differential realization 



Jo = -i4r, J± = e±i4 



(3.2) 



where the auxiliary variable <fi G [0, 2tt) facilitates their closure and the two functions 
F(x), G(x) G C are subjected to constraints of the form 

d ^ = \-F\ ^± = -FG, xER. (3.3) 
ax ax 

Note that we have here J_ ^ thereby inducing an sl(2,C) algebra rather than 
so(2,l), which is consistent with J_ = 



In the case of either sl(2,C) or so(2,l), the irreducible representations are fur- 
nished by 

J \km) = m\km), m — k, k + 1, k + 2, . . . , 

J 2 \km) = k(k-l)\km), (3.4) 

which are essentially of the type . The Casimir operator J 2 corresponds to 

J 2 = Jq =f Jo — J±J T 

- &-{w + l) F ' + 2 % G '- G2 -l- (3 - 5) 

Looking for representations that are 

\km) = ^km(x, 0) = i) km (x)^==, (3.6) 

v 2ir 



where k > and m — k + n, n — 0, 1, 2, . . . , we find ipkm{x) to satisfy the 
Schrodinger equation 

- Vim + Kn^fcm = ~ (k ~ lf) km . (3.7) 

In Q3.7| ), the one-parameter family of potentials V™ is given by 

V m = - (m - 1) (m + 1) F' + 2mG" + G 2 

= - (m - |) (m + |) (l - F 2 ) - 2mFG + G 2 , (3.8) 

where Eq. ( |3.3|) has been used. These potentials share the same real energy eigen- 
values 

E^ = -(m-n-\)\ (3.9) 

thus producing a potential algebra. 

Equations ( |3.3j ) can be solved for the functions F and G to get a quite complete 
realization of the sl(2,C) algebra. The results obtained by us may be summarized 
as follows: 

I : F(x) = tanh(x — c — i7), G(x) = bsech(x — c — i'j), 

II: F(x) = coth(x — c — i7), G(x) = 6cosech(x — c — i7), (3.10) 

III : F(x) = ±1, G(x) = be^ x , 
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where b — &r + ibj, bj G R, and — j < 7 < f . These lead to potentials 



X rn 



(b R + ibi 



m + 



sech 2 (x — c — i'y) 



II : V. 



Ill 



— 2m{b R + ibj) sech(x — c—'vy) tanh(x — c — ij), (3.11) 
(b R + ibi) 2 + m 2 — j cosech 2 (x — c — i'j) 

— 2m(b R + ibj) cosech(a; — c—'vy) coth(x — c — 17), (3.12) 
V m = (b R + ib I ) 2 e T2x T2m(b R + ib I )e Tx . (3.13) 



In this way the group theoretical approach of the potential algebras can be ex- 
tended to non- Hermit ian Hamiltonians (a subclass of which forms the PT-symmetric 
ones) by a simple complexification of the real algebras considered for Hermitian 
Hamiltonians. 

4 sl(2,C) Potentials in SUSYQM 



To realize sl(2,C) potentials from the supersymmetry-inspired Eqs. ( |2.20|) - (|2.23|) , 
we notice that V m can be considered as a special case of the complex potential 
\/(+) = W 2 — W' + E given by (|2.8| ), corresponding to the choice of the complex 
superpotential 

W = (m - I) F - G, 



and the real energy 

E = E R = -(m-\) 2 . 
Inserting fl4.1|) and ([4.2|) into the definition of , we indeed get 
V^=V m = [(m - \) F - G] 2 - [(w- I) (l-F 2 ) +FG 

= ~ ( m " \) ( m + \) i 1 - p2 ) - 2mFG + ° 2 > 
which coincides with the sl(2,C) form ( |3.8| ). 

The potential ' = W 2 + W + E of the superpartner is now 



(4.1) 
(4.2) 



— m — ~ 



(4.3) 



(m - I) F - G] + [(m - I) (l - F 2 ) + FG - 
-(m-|) (m- I) (l - F 2 ) - 2(m - 1)FG + G 2 

Kn-i. 



(4.4) 



From ( |3.9| ), it is obvious that we are here in the case where if_ has one level less 
than H + and Eq. (|2.9|) applies. 

In Table 1 we have displayed the various forms of the complex superpotential W 
for different solutions of F and G summarized in (|3.10|) . 

It is interesting to discuss the results corresponding to the choice 7 = 0. While 
for bi = 0, the superpotential along with the partner potentials reduce to their real 
forms which is only expected, the possibility 5r = is worth taking a look. For 
case I, W simply becomes 

W = (m — tanh(x — c) — ibi sech(x — c), (4.5) 

leading to 

y(+) = y m = (^—b 2 j — m 2 + |^ sech 2 (x — c) — 2imfe/ sech(x — c) tanh(x — c). (4.6) 
Its superpartner can be read off readily from (|4.4|) and is 



' = V m _i = —bj — (m — l) 2 + j sech 2 (x— c)— 2i(m— 1)6/ sech(x— c) tanh(x— c). 

(4.7) 

Note that both and turn out to be PT symmetric. For completeness we 
give the solutions for /, g, Vr , and Vj . These are 

/ = (m- i)tanh(ar-c), (4.8) 

g = — bi sech(x — c), (4.9) 

V^ +) = (-6 2 -m 2 + i)sech 2 (x-c), (4.10) 

Vj +S} = — 2mbi sech(x — c) tanh(x — c), (4-H) 

-^-(m-l) 2 + ilsech 2 (x-c), (4.12) 



V R 



V} ] = -2(m- l)6/sech(x-c)tanh(x-c). (4.13) 

A particular case of the above scheme corresponding to m = 1 was derived 
by Bagchi and Roychoudhury |0J, who showed that the PT-symmetric combination 
of (|4.10|) and (|4.11|) has energy levels that are negative semi-definite and, except for 



the zero-energy state, coincide with those of the sech 2 potential resulting from (|4.12[) 
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Similarly we can deal with cases II and III, both of which are PT non-symmetric, 
for the choice of parameters 7 = and 6r = 0. While case II gives the coth, cosech 
version of (|4.8|) -( |4TT3|) , case III leads to the complexified Morse potential. The 
results for W, f, and g are 
Case II for 7 = 0, b R = 0: 

W = (m — I J coth(a; — c) — ibj cosech(x — c), 
/ = (m — |) coth(x — c), 

g = —bi cosech (x — c), (4-14) 
Case III for 7 = 0, b R = 0: 



w = 


± (m — 


\) 


f = 


± (m — 


i) 


9 = 


-b T e* x . 





(4.15) 

5 A PT- Symmetric Potential in Terms of Weier- 
strass p Function 

In this section we make a specific attempt to search for a PT-symmetric potential 
described by Weierstrass p function. In Ref. 0, Andrianov et al. took V R = 0, 
Vj = to analyze complex transparent potentials belonging to the set ( |2.20| )- 
( |2.23| ). Here we consider an equally viable possibility by setting V R = 0, Vj = 0. 
This case is indeed nontrivial since other possibilities are either related to Andrianov 
et al. conjecture or the present one under f —*■ —f and g — ► —g. 



While Vj i+) = results in 



/ = f, (5-1) 



leading to 



r(+) 3g^_f 
4# 2 2g 



Vk + > = -^-WZ-9 2 + E R , (5.2) 



- -$ + C^ 2 + ^ (5-3) 
Vt ] = 2gf, (5.4) 
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V R ^ = gives us the solution for g in terms of the differential equation 

-. = = = ±dx, (5.5) 

Ig (|(?3 - AE R g + a) 

where a represents a constant of integration, which we take to be non-zero. 

Writing y(g) = g (|g 3 — AE R g + aj in the form y(g) = a g 4 + 4ai# 3 + 6a 2 g 2 + 
4a 3 g + a 4 , we have a — |, a\ — 0, a 2 = — \Er, a 3 = |a, a 4 = 0. We next define 
quantities gi and g 3 as 

g 2 = a a4 — 4a!03 + 3a 2 = ^E R , (5-6) 

#3 = ao°204 + 2aia 2 a3 — a 3 , — 0003 — a 2 ^ = J^-E^ — ^Q 2 - (5-7) 

Let 

« = T fe/wr 1/2 ^ ( 5 - 8 ) 

"'go 

where g is any root of the equation y(g) = 0. We identify go as go = 0. 
Let us substitute £ = ^ and # = | to rewrite (|5.8| ) as 

2 = y (4a 3 r 3 + 6a 2 r 2 + a ) dr. (5.9) 

The second term in the integrand can be removed by effecting the transformations 

U*+lE R ), (5.10) 



CI3 a \ 3 

03 a V 3 / 
As a result, z turns out to be given by 

/°° / \ —1/2 

(4a 3 -0 2 <7-0 3 ) ^- (5-12) 

Now from the theory of elliptic functions [[TIJ, we can read off 

s = p(z;g 2 ,g 3 ), (5.13) 



where p(z; g 2 , #3) is Weierstrass elliptic function with g 2 and #3 as invariants. Equa- 
tion (|5.11|) yields 

. (5.14) 



9= A 



p(z;g 2 ,g 3 ) + ^E R 
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Hence we deduce 



Vi 



(-) 



(+) 



P(z) + \E R 
1 



2 ' 



:E f 
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p(z) + \E R 



2 ■ 



(5.15) 
(5.16) 



In deriving Q5.16), use has been made of the differential equation satisfied by p{z) 



and of its consequence 



p'\z) = Ap i (z) -g 2 p(z) -g 3 , 



p"{z) = 6p 2 ( Z ) - §0 2 . 



(5.17) 



(5.18) 



We notice that while V i+) = is pure real, = iVj { ~ } is a pure imaginary 
potential. 

Comparing fl5.5| ) with ( |5.8| ), we see that z = ±x + c, where c is an integration 
constant. Since p(z) and p'(z) are respectively even and odd functions of z, the two 
solutions obtained by taking either sign correspond to each other by a mere change 
of signs of the integration constants a and c. Thus without any loss of generality we 
can take z = x + c only. 

We may distinguish between the non-degenerate and degenerate cases of Weier- 
strass p function [|K|. 

In the non-degenerate case, the roots of the cubic equation 



4cr - g 2 (? - g 3 = 
are all distinct and the corresponding discriminant 

,2 



(5.19) 



D = gi- 27g : 



a 
48 



6AE R - 9a' 



(5.20) 



is non-vanishing. It is either positive or negative according to whether \a\ < ^E^ 2 
or \a\ > ^Eft 2 , the former case occurring only for En > 0. 

By using numerical studies, we showed that in the D < case, wherein the 
Weierstrass function has a pair of complex conjugate primitive periods 2uj, 2oj' = 
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2u*, V { R +) and v} } go to -oo for some z values because the denominators in ( |5.15| ) 
and Q5.16 ) vanish for such values. On the contrary, in the D > case, wherein 
the Weierstrass function has a pair of primitive periods 2u, 2a/ with uo real and u' 
imaginary, we obtain well-behaved potentials defined on the interval < z < 2uo or 
— c < x < 2u — c . The potential is a single-well potential, singular at z — > 
and z — > 2a; (where it behaves as I/2 2 and l/(z — 2u) 2 , respectively), and symmetric 
around its minimum at z — u, whereas Vj~^ vanishes at z —>■ and z — > 2uo, and is 
antisymmetric around z — ou. Hence, the potential y(~) = iVj "* is PT symmetric 
provided the parity operation is defined with respect to a mirror placed at z = u or 

X = UJ — c. 

In Fig. 1, the functions and Vj~^ are displayed in terms of z for E R = \/3 and 
a = 4y / 2/y / 3, corresponding to the invariants g 2 = 4, g 3 = 0. For such values, the 
cubic equation ( |5.19| ) has the three real roots e\ = 1, e 2 = 0, = — 1, and the real 
primitive period is given by 2uj = 2x* where x* = a/tt r(5/4)/r(3/4) ~ 1.311 fTT| . 
The minimum of V R + ^ is equal to 6 (l — . 

In the degenerate case, at least two of the roots of ( |5.19| ) are equal, meaning 
that D = 0. This condition imposes that E R > 0, a = ±^E R 2 , g 2 = ^E R , and 
g 3 = — -^=E R . We are then in a case where the real period becomes infinite and p(z) 
reduces to [jlCf 



In consequence we have 



E R \ + cosech 2 ( J E R z 



fj 



1 + § cosech 2 U/E^z) 



(5.21) 



(5.22) 



from which we obtain 



(+) 



R 



4 I 3 

-Er < 1 H — cosech 2 
3 1 2 



Erz) - 



I cosech 2 



(-) 



±6£, 



cosech 2 (s/Er, z^j coth [\/E R z^j 



1 + § cosech 2 (v^fl^) 



2 ( ' 



(5.23) 



(5.24) 
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defined on the interval < z < oo or — c < x < oo. As such \ given by 
= iV^, is a non-PT-symmetric potential. We also note that as z — > 0, 
V^ +) ~ 2/z 2 while for z — > oo, V^ +) ~ 24£ fl exp(-2 % /i^). These are reasonable 
boundary conditions, the behaviour of proving to be repulsive. 

An example is displayed in Fig. 2 for E R = y/3 and a = 8/3 1 / 4 , corresponding to 
the invariants g% = 4 and g% = — 8/3 3 / 2 . For such values, the cubic equation ( 5.1 9| ) 
has the three real roots e% = e 2 = l/VS and = —2/y/3. 

As discussed in Sec. 2, the spectra of the supersymmetric partners H + and H_ are 
alike except for the ground state. In the non-degenerate case, whenever E R > and 
\a\ < |-E^ 2 , the Hermitian Hamiltonian H + has an infinite number of (unknown) 
discrete positive-energy levels. The same is true for the PT-symmetric Hamiltonian 
but in addition the latter has a normalizable eigenfunction ipQ ^ corresponding 
to E = 0. From (|2~T0| ), (gig ), and fl5l]), it is given by 



(5.25) 



By using ( |5.14| ), its modulus can be written as 

\K\ 



p(z;g 2 ,g3) + -E R 



-1/2 



(5.26) 



Hence it vanishes at z = and z = 2u, and is regular in between, showing that 
ipQ \x) is indeed normalizable on (— c, 2uo — c). 

Similarly, in the degenerate case, i.e., for Er > and \a\ — ^E R 2 , the Hamilto- 
nians H + and H_ have both a continuous spectrum of unbounded positive-energy 
states. In addition, H_ has an unbound zero-energy state, whose wave function is 
still given by (|5.25 ). From (|5.14j) and (|5.21|) , we indeed obtain 

-1/2 



cosech z 



E R z) + - 



(5.27) 



which vanishes for z = and goes to ^ \/Hf^ f° r 



oo. 



As a final point, it is worth noting that the known zero-energy eigenfunction of 
H is an eigenfunction of the Hamiltonian — -pr + V"(~)(a:) with energy .Er. 
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6 Conclusion 



In the present paper, we both constructed some new PT-preserving or non-PT- 
preserving complex potentials and analyzed some known ones from a SUSYQM 
viewpoint extended to deal with non-Hermitian Hamiltonians. 

To start with, we considered three sets of complex potentials, recently derived 
by a potential algebraic approach based on the complex Lie algebra sl(2,C) |J, and 
proved that they can be generated as well from a complex superpotential and a 
pair of supercharges that are not related by Hermitian conjugation. This extends to 
the complex domain the well-known relationship between SUSYQM and potential 
algebras for Hermitian Hamiltonians, resulting from their common link (12| with the 
factorization method il3fl and Darboux transformations HT4 . 



In a second step, we analyzed the special case of the extended SUSYQM the- 
ory H wherein the starting potential is real and its partner imaginary. This allowed 
us to build for the first time a pair of complex partner potentials defined in terms 
of Weierstrass elliptic function. The PT-symmetric imaginary partner potential has 
one known eigenvalue equal to the factorization energy Er and corresponding to a 
bound state. When the Weierstrass function degenerates to a hyperbolic one, the 
imaginary partner potential becomes PT non-symmetric and its eigenvalue Er cor- 
responds to an unbound state. We have therefore constructed two new quasiexactly 
solvable complex potentials 0, 0]. 
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Table 1: The functions F, G, and the superpotential W corresponding to the three 
cases considered in the text. 



Cases F G W 

I tanh(x — c — i7) b sech(x — c — ij) (m — |) tanh(a; — c—'vy) 

— 5sech(:r — c — vy) 

II coth(x — c — 'vy) b cosech(:r — c — 'vy) (m — |) coth(x — c — 'vy) 

— 6cosech(a; — c — i^y) 

III ±1 be TX ± (m - |) - bc~^ x 
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Figure captions 

Fig. 1. (a) and (b) Vj~^ in terms of z = x + c for the non-degenerate case of 
Weierstrass p function, E R = y/3, and a = 4^2/ \/3- 

Fig. 2. (a) and (b) in terms of z = x + c for the degenerate case of 

Weierstrass p function and E R = a/3. 
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Figure 2 
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